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1. Introduction 

In recent years there has been an increasing interest in the geometric treatment of non- 
holonomic mechanical systems, see, e.g., [2,3,5,6, 10, 12, 19,28,29, 31,38-40]. In particular, 
it has been recognised that the Hamiltonian formulation of such systems can be stated in terms 
of an almost-Poisson bracket, that is, a biderivation of functions of phase space, antisymmetric 
in its arguments but which does not necessarily fulfil the Jacobi identity (see, e.g., [1, 1 1, 36]). 
Therefore, for researchers in this field, it seems to be usual the conceptual association of 
the Hamiltonian formulation of non-holonomic mechanical systems with almost-Poisson 
structures. 

On the other hand, there exist non-holonomic systems which, after certain reductions 
are performed, admit a Hamiltonian formulation after a "rescaling of time" is carried out, by 
means of rescaling factors (sometimes called invariant measures) of the reduced vector field 
of the system. This is the case for the so-called LR systems, which are systems formulated on 
compact Lie groups endowed with a left-invariant metric and right-invariant non-holonomic 
constraints. After a rescaling of time, their corresponding reduced systems become integrable 
Hamiltonian systems describing geodesic flows on unit spheres [24]. In [9], a necessary and 
sufficient condition for the existence of an invariant measure for the reduced dynamics of 
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generalized Chaplygin systems of mechanical type is given. Another recent work on this 
line is [43]. For a classic treatment of the theory of Chaplygin's reducing multiplier, see 
Section 111-12 of [33]. Thus, it could be conceptually associated as well the existence of 
specific rescaling factors for these reduced systems with the possibility of formulating them 
in a Hamiltonian way. 

In addition, Borisov, Mamaev and Kilin [7,8] have recently found a Poisson structure for 
each studied case of reduced non-holonomic systems, such that the reduced system becomes 
Hamiltonian, with respect to such a structure, after a rescaling, the Hamiltonian function being 
the reduced energy. The examples treated by them are classical in the literature, consisting 
mainly of rolling bodies without slipping, namely a rigid body of revolution rolling on a 
plane, in particular the Routh's sphere (see Section 4.2), the rolling disk (to be treated in 
Section 4.1), the motion of a homogeneous ball on a surface of revolution (called sometimes 
Routh's problem, see, e.g., [42,43]), and other cases. There is a strong emphasis in these 
references in the sense that the Poisson structure for each case can be found after a rescaling 
of time of the reduced vector field. 

Our primary motivation for this work was to understand the origin of the two integrals of 
motion appearing in the mentioned problem of a ball rolling without slipping inside a surface 
of revolution, which are not given, in general, in an explicit form but being related to the 
solutions of a system of first order non-autonomous differential equations [25,35,42]. This 
also happens in the other mentioned cases. The results of [7, 8] suggest that such systems can 
be interpreted as the equations providing a set of functionally independent Casimir functions 
of the Poisson structure they find for each specific case. Therefore, it seemed to be worth 
investigating further such Poisson structures, in particular to clarify their domain of definition 
and basic properties. Let us note that another recent approach, devoted to the study of Poisson 
structures which can be associated to never vanishing vector fields on manifolds of arbitrary 
dimension d > 2, with fibrating periodic flows, is given in [23]. 

It follows that the previously mentioned Poisson structures have a rather peculiar form. 
In particular, the associated characteristic distributions have rank two in the open sets of 
the reduced spaces considered in [7,8]. This property implies that such Poisson structures, 
when multiplied by a never vanishing function, are again Poisson structures of the same type. 
The immediate consequence is that the above mentioned reduced non-holonomic systems 
are already Hamiltonian with respect to one of these Poisson structures without any need of 
rescaling. 

Other interesting result is that, in the cases studied, the Poisson structures obtained can 
be extended from their original domains of definition, namely (open sets of) semialgebraic 
subvarieties of M 5 , to an open set of the ambient space. Such extended Poisson structures 
become zero only at the so-called singular equilibria of the reduced systems. Moreover, the 
existence of these (extended) Poison structures, from an algebraic point of view, is only caused 
by the existence of integrals of motion of the reduced vector field related to the solutions of 
the mentioned systems of first order differential equations. 

This paper is organized as follows. In Section 2 we briefly review some notions of 
Poisson geometry and in particular, of Poisson structures of rank two. In Section 2 we show 
the explicit expressions of certain bivectors in K 4 and R 5 , determined up to a non-vanishing 
factor function, by choosing the 1 -forms in their kernels to have a specific form, and we 
prove that they are in fact Poisson bivectors of rank two. Section 4 is devoted to show the 
application of the previous results in specific examples, namely, the rolling disk, the Routh's 
sphere, and the ball rolling on a surface of revolution. We will use the formulation of [ 1 8] , [ 1 7] 
and [25], respectively, of these problems, rather than that of [7, 8]. However, we point out the 
equivalence of both treatments in the last case. We also treat the special case of a ball rolling 
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inside a cylinder. Finally, we end with some conclusions and an outlook for further research. 
2. On Poisson structures of rank two 

For the sake of completeness and in order to fix some notations, we will recall some well- 
known notions on Poisson manifolds, and in particular, we will focus on Poisson structures of 
rank two. For more details see, e.g., [26]. 

Given a differentiable manifold M, a Poisson structure on M is defined by an 
antisymmetric bilinear map {-, ■} which is a derivation on both of its arguments, satisfying 
moreover the Jacobi identity. A manifold M endowed with a Poisson structure is called a 
Poisson manifold. 

Thus, it is possible to associate to each function / a unique vector field Xf such that, for 
any other function g, we have Xf-g = {f,g}. The vector field Xf is called Hamiltonian vector 
field associated to the Hamiltonian function f. This association defines an homomorphism 
of the Lie algebra (C°°(M), {•,•}) onto the Lie algebra of vector fields in M. A Casimir 
function or Casimir for short, is a function c such that X c = 0. The Poisson structure is called 
non degenerate if only the constant functions are Casimir functions. 

Moreover, on every Poisson manifold, there exists a unique twice contravariant antisym- 
metric tensor field (called bivector field for short) A such that {/, g} = A(df, dg) for every 
pair of functions (/, g). This tensor field is called the Poisson tensor of the structure, and the 
manifold M, endowed with its Poisson structure, will be denoted (M, A). The existence of 
such a tensor field is due only to the antisymmetry and derivation properties of the Poisson 
bracket. The fulfillment of the Jacobi identity for the Poisson bracket is equivalent [27] to 
the vanishing of the Schouten-Nijenhuis bracket of A with itself, [A, A] = 0. The Schouten- 
Nijenhuis bracket [34, 37] is the unique extension of the Lie bracket of vector fields to the 
exterior algebra of multivector fields. Some of its properties are 

[P,Q] = -(-i)< p - 1 ><*- 1 >[Q,P] 

[P, Q A R] = [P, Q] A R + {-l) ip - 1)q Q A [P, R] (1) 
[P A R,Q] = P A [R, Q] + (-l)( q -V r [P, Q] A R 

where P,Q,R are completely antisymmetric contravariant tensors of degree p,q,r, res- 
pectively. For more details and properties on the Schouten-Nijenhuis bracket see, e.g., 
[13, 30, 34, 37] and references therein. 

Take a local chart of M, with domain U and associated local coordinates (x\, . . . , x n ), 
where n = dimM. We will denote by Ay, (1 < i,j < n) the components of the Poisson 
tensor A in the previous chart. The expression of the Poisson bracket of the restriction of the 
two functions /, g to U, also denoted by /, g, reads 

df df 



if, 9} = A y 



dxi dxj 



where summation in the repeated indices is understood. In particular we have {xi , Xj } = A y . 
The Poisson tensor admits the local expression 

in these coordinates. 

Given a Poisson manifold (M, A), it can be defined the fibered morphism A' : T*M — > 
TM such that for any pair of 1-forms a, (3, (A^a),/?) = A(a,/3). The image of the 
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morphism A', C = A"(T*M), is called the characteristic distribution of the Poisson 
structure, and the characteristic space on x G M is the vectorial subspace C x — A X (T*M) 
of T^Af. The ranfe of the structure on the point x is the rank of A x , i.e., the dimension 
of C x . Note that the annihilator of the characteristic distribution, i.e., C° = {[3 G 
A X (M) | A(/?,a) = 0,Va£ A 1 (M)}, is kcr A J , and we have rankA| + dimker Ajj. = n, 
for all x e M. In general, the rank of the structure varies with a; and thus C is not in general 
asubbundle of TM. 

Consider now a Poisson manifold (M, A), dimM = n, such that in the domain of a 
local chart (U, <f>) the structure has constant rank equal to two. The Theorem 1 1.5 of Chapter 
III in [26] (or Corollary 2.3. in [41]) assures us that the associated local coordinates, denoted 
(x, y, zi, ... , Zn-2), can be chosen such that for 1 < k , I < n — 2, 

{y,x} = l, {x,z k } = 0, {y,z k } = 0, {z k ,zi} = Q. (3) 

We are now in a position to prove a simple result, but important for our purposes here: 

Proposition 1 Let (M, A) be a Poisson manifold of (locally) constant rank equal to two. 
Then, for each never vanishing smooth function a <G C°°(M), (M, a A) is a Poisson manifold 
of (locally) constant rank equal to two, with the same characteristic distribution. 

Proof 

We have to prove that the Schouten-Nijenhuis bracket [a A, aA] vanishes, the other 
needed properties being obvious. From the paragraph 18.8 of Chapter V of [26], we have 
that 

[aA, aA] = 2aA J (da) A A. 

It suffices to compute the previous expression on a coordinate neighbourhood like that 
described in the previous paragraph, with respect to the Poisson tensor A [30]. We have 

[aA, aA](dx, dy,dz k ) = 2a(A tt (da) A A)(dx, dy, dz k ) = 0, 1 < k < n - 2 , 

because z k are Casimir functions of A, and dz k enter at least once as argument of A in all 
terms of the previous expression. For other possible arguments the expression vanishes by the 
same reason. 

Example 1 Let M be a n-dimensional manifold and X, Y two vector fields such that for 
all x <G M, the Lie bracket [X, Y] x belongs to the subspace of T X M generated by X x and 
Y x . Then, A = X A Y is a Poisson tensor of rank two except where X and Y are linearly 
dependent. This is easily seen by deducing from the properties of the Schouten-Nijenhuis 
bracket (1) the relation [X A Y, X A Y] = 2X A Y A [X, Y], see also [1, 13]. 

Remark Note that it is essential in Proposition 1 the assumption that the initial bivector is 
Poisson, which assures the existence of local coordinates satisfying (3). The existence of a 
bivector whose rank is always two is not enough to conclude that it is a Poisson bivector. A 
simple counter-example is the following. Take M = M. 3 , with coordinates (x, y, z). Let X, 
Y be vector fields given by 

X- — - — Y-— + X— 

dx ^ dz ' dy X dz 

Then, A = X AY = J^A^Msan everywhere rank two bivector but is not Poisson, since 
[X,F]-2^and 

[ A,A ] = [XAY,XAY ] =^ x A^A§-. 
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The vector fields X, Y and [X, Y] in this example close on a Lie algebra isomorphic to the 
Heisenberg-Weyl Lie algebra h(3), see, e.g., [16]. 



3. Some Poisson structures of rank two in R 4 and R 5 



We will construct in this Section some Poisson structures of rank two in R 4 and R 5 
by imposing that the kernel of the corresponding bivectors consists of a set of two and 
three specific 1 -forms, respectively. Such 1 -forms will determine codistributions which are 
integrable in the sense of Frobenius. We will prove that the resulting bivectors are in fact 
Poisson. 



3.1. Some Poisson structures of rank two in R 4 

Consider the Euclidean space M 4 , with coordinates (x\,x 2 , x 3 , £4). The equations of motion 
of the reduced non-holonomic systems encountered in the examples are observed to have 
integrals of motion which are related to the solutions of a system of differential equations of 
the type 

^ = h 3 (xi, x 3 ,x 4 ) , ^ = hi{xi,x 3 ,Xi) , (4) 
dxi dxi 

where h 3 , /14, are two given (smooth) functions of their arguments, which do not include X2- 

We consider the system (4) as the Pfaffian system '8± = 0, 82 = 0', where the 1-forms 61, 62 

in R 4 , are given by 

6\ = — hs(x\, xs,X4)dx\ + dx 3 , 9 2 = —h4(xi,x 3 ,X4)dxi+dx4. (5) 

These two 1-forms determine a codistribution integrable in the sense of Frobenius [26], since 
there exist a set of four 1-forms Aj such that d9i = A\ A 6j for i, j — 1,2. For example, we 
can take 

A[ = -^—dx! , Af = -^—dxi , A2 = ^— dxi , A^ = -^—dxi , (6) 

axz 0x4 0x3 0x4 

in order to satisfy the integrability condition. Thus, there will exist (locally) functions c\ , c 2 
such that 6i = dc i7 i = 1,2. The subvarieties solution of the Pfaffian system '81 — 0, 82 = 0' 
are defined by the equations Cj = hi, where bi are constants, i = 1,2. 

More specifically, in the actual examples, the system (4) takes the form of a non- 
autonomous first order system of linear differential equations 

^1 = an(xi)x 3 + a 12 (xi)x4 , ^— ^ = a 2 i(xi)x 3 + a 22 (a;i)x4 , 
dxi dxi 

or, written in matrix form, 

±( X " )=A( X1 )( X3 V (7) 
dxi \ x 4 J \ x 4 J 

where 



an(xi) ai 2 (xi) 
021(2:1) a 2 2(xi) 

The previous functions Cj can be identified with the initial conditions of the solution of (7). In 
fact, such a solution can be expressed as x = g(x\)c, where x = (x 3 , x 4 ) T , c = (ci, C2) T , 
and g(x{) is a GL(2, R)-valued curve (SL(2, R)-valued curve if tr A(xi) = for all x\), 
solution of the right-invariant matrix system (see, e.g., [14, 15]) 

^-g- x =A{x{). (8) 
dx\ 
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Then, c = g~ 1 (xi)x gives the desired functions: with a slight abuse of notation, we have 

dc = (dg~ 1 )x + g~ 1 dx. = —g~ 1 dgg~ 1 x + g^ 1 Axdxi = —g~ 1 dgg~ 1 x + g~ x dgg~ x x = 0, 

where we have used that dg^ 1 = —g~ 1 dgg~ 1 and Adxi = dgg^ 1 . However, note that the 
solution of (8) cannot be expressed in an explicit way in the general case, and therefore, the 
functions c\, c 2 cannot be explicitly written in general. 

Now, we impose that the 1 -forms (5) generate the kernel of the bivector in R 4 

A = E A ^ A 7^- ^ 

The resulting bivectors will clearly have rank two. Moreover, they are Poisson, according to 
the following result 

Theorem 1 Consider in R 4 a bivector of type (9), such that A^i) = 0, A^(9 2 ) = 0, where 
Q\, 9 2 are given by (5). Then the bivector is of the form 



A= -Ki 2 U AV , 



where 



d d d d 

U=j—, V= — + h 3 — + h 4 —, (10) 
0x2 oxi 0x3 0x4 

and A12 € C°°(M 4 ). Each of these bivectors is Poisson, and of rank two on points where 
A12 + 0. 

Proof 

The case of A12 = is trivial. We will assume A12 ^ in the domain of interest. Take 
A and 6\, 6 2 , as stated. The conditions A"(#i) = 0, A" (#2) = give rise to an algebraic 
system for the six independent functions Ay, which can be easily solved for five of them, in 
terms of the remaining one and the functions entering into the 1 -forms. We choose A 12 to be 
the undetermined function. Then the solution reads 

A13 = A w = A 34 = , A 23 = -A 12 /i 3 , A 24 = -A 12 /i 4 , 

thus the resulting bivectors are as claimed. To see that each of them is Poisson, consider the 
bivector of the family with A12 = — 1, i.e., U A V. This bivector is of the form given in 
Example 1, and [U, V] = 0, thus U A V is Poisson. It is moreover of rank two, therefore by 
Proposition 1, the claim follows. 

Remark Note that the vector fields U, V of the previous Theorem satisfy 9 l (U) = 6 t (V) = 0, 
i = 1,2, which in principle might seem a stronger condition than that the bivector (9) 
annihilates the 1 -forms 61,62- 

Now, given a (Hamiltonian) function H <E C°°(M 4 ), the Hamiltonian vector field Xh 
with respect to a Poisson structure of the family described on Theorem 1 takes the form 

X H = K\dH) = A 12 [(VH)U - (UH)V] , (11) 

where U and V are given by (10). Obviously, H is a first integral of Xh, since XhH = 
A(dH, dH) = 0. Other two first integrals are the functions Cj such that dci = 6i, since by 
construction X H (ci) = A(dH, dc{) = A(dH, 6i) = 0, i = 1, 2. These two first integrals are 
common to all Hamiltonian vector fields of type (11). 

On the other hand, given a specific vector field X in M 4 , which is recognized to be of 
the form (1 1), it could be regarded as a Hamiltonian vector field with respect to one specific 
Poisson structure of the family described in Theorem 1 . 



7 



3.2. Some Poisson structures of rank two in R 5 

We will treat in this Section analogous questions to that of the previous Section, but now in 
the Euclidean space R 5 , with coordinates (a;i, X2, £3, X4, £5). 

The motivation is that typically, the reduced orbit spaces for the non-holonomic problems 
of interest, are semialgebraic varieties of R 5 , essentially determined by the zero level set 
of a function cfr € C°°(R 5 ), quadratic in its arguments, which are moreover subject to 
certain constraints. More specifically, in the examples it will have the form 4>(x) = 0, with 
(f){x) = x\-\-x\ — (1 — x\ )x 5 , \xi \ < 1, and x 5 > 0, or with (f>(x) = x\ + x\ — 4xix 5 , X\ > 0, 
and x 5 > 0. However, for what follows <p can be in principle any differentiable function in 
R 5 . 

We will consider then the Pfaffian system '0 O = 0, 9i = 0, 62 = 0', where 6 = d(f> and 
9i, 02 are 1 -forms in R 5 whose coordinate expression is again (5). These three 1 -forms also 
determine determine a codistribution integrable in the sense of Frobenius in R 5 , because we 
have again d6i = A? A 6j with (6), i, j = 1, 2, and d9 = d 2 (f> = 0. 

We impose now that kcr A' = span{^ , 61, O2}, where A is the bivector in (some open 
set of) R 5 

A - £ A ^ A ^- (12) 

The resulting bivectors are again generically of rank two and Poisson, as follows 

Theorem 2 Consider in R 5 a bivector of type (12), such that A*(0 O ) = 0, A*(#i) = and 
A" (O2) = 0, where 9o = d(j), and Q\, 9 2 are given by (5). Then the bivector is of the form 

A = f[(Z<f>)U AV + Y AZ], (13) 

where 

d d d d d 

0x2 0x1 0x3 0x4 0x5 

Y = (U(p)V - (V(f>)U , (15) 

and f € C°° (R 5 ). Each of these bivectors is Poisson, and of rank two on points where / ^ 0. 



Proof 

Once more, the case of / = is trivial, thus we will assume again that / ^ in the 
domain of interest. Take A, 6q, Q\ and 62 as stated. The idea of the proof is similar to that of 
Theorem 1 . First of all, since the kernel of A" has generically dimension three, then the rank 
of A" is two. The conditions A"(0o) = 0, A^fli) = and A i (9 2 ) = give rise again to an 
algebraic system for the functions Ay , out of which all can be solved for except one of them, 
namely A12, which we will write as — (dtfi/dx^f. The solution then reads 

A13 = A M = A 34 = , A 23 = fh 3 g^ , A 24 = fhi-Q^ . 

A15 = Jt; — j A 35 = j/13- — , 
0x2 0x2 

a n dfi A n f 9(f> . d<f> . d<f> 

A45 = fhAj— , A 25 = -/ — + h 3 — + h A — 

OX2 \OXi OX3 0x4 

thus the resulting bivectors take the stated form. To see that each of them is Poisson, consider 
the bivector of the family with / = 1, i.e., A = UAV + YAZ, where U = (Z<f>)U. We 
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have to show that the Schouten-Nijenhuis bracket of A with itself vanish, i.e., [A , A ] = 0. 
By linearity and using the first property of (1) we have 

[Ao, A ] = [F A V,U A V] + 2 [17 A V, Y A Z] + [Y A Z, Y A Z] 

By Example 1 we know that \U A V,U A V] = 2U A V A \U, V] and analogously, 
[Y A Z, Y A Z] = 2 Y A Z A [Y, Z]. Now, using again the second and third properties of 
(1) we can write 

[U A V, Y A Z] = V A Z A [U, Y] - U A Z A [V, Y] 

+ Y A V A [77, Z] - Y A U A [V, Z] 
We have to calculate now some Lie brackets. We have [U, V] — [V, Z] — [U, Z] = but 

\u,v] = -[v(Z4>)]u, [y,z] = -[z{U4>)]v + [z(Y<f>)]u 

[u,y] = {Z4>)[u{u4>)]v - {{z<}>)[u{V4>)] + (U(f>)[v{z<i>)] - {v<f>)[u{z<p)]}u 

[V,Y] = [V{U4>)]V- [V{V<j>)]U, [U,Z] = -[Z(Z^)]U 

Then, summing up, we have 

[A , A ] = 2 U A V A Z{{Z<f>){\V, U]4>) + (U^Z, V\<j>) + {V<t>)(p, Z]4>)} = . 

Since the rank of any of the A, and in particular A , is two, applying Proposition 1 ends the 
proof. 

Remark Note that the vector fields U, V, Y and Z of Theorem 2 satisfy 0i(U) = 9i{V) = 
6i{Y) = 9i{Z) =0,i = 1,2, 9 (Y) = Y{4>) = and (UAV)<p-Y = 0. These requirements 
might seem a priori to be stronger conditions to that imposed in the Theorem. 

If we are given now a (Hamiltonian) function H e C°°(M 5 ), the Hamiltonian vector 
field Xh with respect to a Poisson structure of the family described in Theorem 2 reads, 
using (15), 

X H = A\dH) = f{[(ZH)(Vcf>) - {Zcj>){VH)] U (16) 

+ [(Z<f>)(UH) - {ZH){U4>)\ V + [{U4>)(VH) - (V<P)(UH)} Z) , 

where U, V and Z are given by (14). By construction H is a first integral of Xh- Other first 
integrals are the functions c» such that dc.i = 9i, as in the previous Section. These two first 
integrals are common to all Hamiltonian vector fields of type (16). 

However, given a specific vector field Xh of type (16), it fixes the specific function / 
and therefore the specific Poisson bivector of the family (13) with respect to which Xh is 
Hamiltonian. 

4. Examples 

In this Section we will show how the preceding results can be directly applied in the cases 
of reduced systems corresponding to specific examples of non-holonomic systems, i.e., the 
rolling disk, the Routh's sphere, the ball rolling on a surface of revolution and its special case 
of a ball rolling inside a cylinder. 

4.1. The rolling disk 

For this example we will follow the treatment and use some of the results of [18], see 
details therein. This problem has been treated as well, e.g., in [4,7, 19,33,35]. Consider 
a homogeneous disk, which rolls without slipping on a horizontal plane under the influence 
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of a vertical gravitational field of strenght g. The resulting non-holonomic system has two 
evident symmetry groups. One is the symmetry group E(2) consisting of translations in the 
horizontal plane and rotations about the vertical axis, and the second is the S 1 symmetry 
consisting of rotations about the principal axis perpendicular to the plane of the disk. 

After these two symmetries have been reduced out, in particular by using invariant theory 
for the reduction of the S 1 symmetry, it is obtained a system giving the evolution on the 
reduced orbit space, which is a semialgebraic variety of R 5 . In particular, the system can be 
restricted to a smooth open subset as it has been done in [18]. 

Thus, consider a reference homogeneous disk of radius r and mass to, lying flat in a fixed 
reference frame with center of mass at the origin. The position of the moving disk is given by 
transforming the position of the reference disk by means of a translation a (e.g., of the center 
of mass) and a rotation A. The tensor of inertia 7 with respect to the principal axes of the 
disk is diagonal, 7 = diag(7i, I\, 7 3 ). Let us call e 3 the vertical unitary vector in the fixed 
frame of reference. We define the unitary vector u with respect that frame as the pre-image 
of — e 3 under the rotation A, u = — A _1 e 3 . The vector s in the fixed disk, rotated by A gives 
the vector in the moving disk pointing from the center of mass to the point of contact of the 
moving disk with the horizontal plane. If we denote u = u — (u, e 3 )e 3 , the relation between 
s and u is s = ru/\u\. We denote by (o>i, w 2 , ^3) the components of the angular velocity 
vector cj of the disk, 

Following [18], after the mentioned symmetry group E{2) is reduced out, the equations 
of motion read 

d(Ioj) 2 dui I du \ . .ds 

— lui x u> — rar — — h to ( ——, s ) s + m(s,u>) — 



dt dt \dt I x ' dt 

+ m(u),s)(u) x s) — rag (u x s) (17) 

du 

— = U X LO 

dt 

which have a first integral given by the total energy of the disk 

77 = i(7w,cj) + i(o> x s,w x s) + rag (s,u) . (18) 

The second of Eqs. (17) expresses the non-holonomic constraint of rolling without slipping, 
i.e., instantaneous velocity of the point of contact equal to zero. 

We recall briefly now how the further reduction of the S 1 symmetry is performed. Let 
us denote by (m, u 2 , u 3 ) the components of u. The S 1 symmetry action consists of rotating 
both vectors u and u simultaneously as mentioned, and it is not a free action since the isotropy 
subgroup of pairs ((0, 0, ±1), (0, 0, w 3 )) is S 1 . Thus, we will use invariant theory in order to 
perform the reduction. A set of invariants for this action is easily constructed [18]: 

<7l = M 3 , <T 2 — U 2 LO\ — U1UJ2 , 03 = U\W\ + U 2 OJ2 , 

<7 4 =u;3, <7 5 = uj\ + col , <r e = ul + ul, (19) 

with the relations 

o\ + <T 3 = (75(76 , CT5 > , (76 > . (20) 

Since u is a unitary vector, we have that <j 6 + af = 1 and |cti| < 1, thus the completely 
reduced orbit space M is the semialgebraic variety of M 5 

M = {(<n,... ,a 5 ) el 5 \<K<t)=0, \<ti\<1, (75 >0}, (21) 

where 4> € C°°(]R 5 ) is the polynomial function (f>{a) = a\ + erf — (1 — a\)a^. However, M 
is not a smooth submanifold of R 5 . The singular points of M are 

n± = {(±1, 0, 0, (7 4 , cr 5 ) e R 5 I (7 4 e R, cr 5 > 0} . (22) 
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The non-smoothness of M is due to the fact that the S 1 action is not free, see [17]. 

The somehow redundant variables (<ti, 02, 03, 04, 05) therefore parametrize the reduced 
orbit space M. The induced system from (17) will be written in terms of the orbit variables: 
simply calculating their time-derivatives, using the equations of motion (17) and that I\ = 
\mr 2 and J 3 = ^mr 2 , we arrive to the following system 

&l = (72 

6 4 crier 2 , r ~ 

°~2 = - (7 3 <7 4 - <7lC7 5 + - g + ACTI W 1 - (7f 

O O 1 (T * 

(73 = -2(72<74 (23) 

2 1 

^ = --3-^(72(73 



1 

2 



(75 = 2ct 2 



A(7i 4 (71C3 4 (73(74 



v /i— ^ 5 (i-(7 2 )2 5i-a 2 y ' 

where A = | fi and the dot means derivative with respect to time. The reduced energy, 



5 r 

obtained from (18), reads 



^T+rHi^+V^- <24) 

Although in principle the expressions (23) and (24) are only defined on M, their 
right hand sides make sense for V = M 5 \({(±1, cr 2 , 03, (74, er 5 ) | ct 2 03 7^ 0} U 
{(ci, o"2 , cr 3 , (74 , (75 ) I |(7i| > 1}), so we will consider this extended domain for the vector 
field X whose integral curves are given by (23) and the reduced energy function E. 

However, if we restrict ourselves to the original domain M, and morever to points with 
|cti I < 1, we can define a smooth open dense subset M C M given by 

{22 ^ 
(<7i, (7 2 , (73,(74,(75) e K 5 I (7 5 = < 1 J ' (25) 

diffeomorphic to R 4 [18]. The induced vector field X and energy E on M can be easily found 
from (23) and (24) by solving for cr 5 . The integral curves of X are the solutions of the system 



(71 = C72 

6 <7l 2 1 a l 
02 = T CT304 - 2 CT 2 _ c 1 

5 1 — af 5 1 - (75 

(73 = — 2(72C4 

2 1 

3 1 — erf 



■crl + Acr^l-cr 2 (26) 



meanwhile 



B ^l4 + ST4 + r* + A ^ <27) 

These expressions are Eqs. (18) and (19) of [18], respectively. 

The reduced vector field X satisfies X(E) = as well as X((f>) = in T>, meanwhile 
X(E) = in M. In addition, X has a family of equilibrium points belonging to the singular 
set n±, called singular equilibria, given by {(±1, 0, 0, CT4, 0) | 04 <G K}, and a family of 
regular equilibria given by the set of constants 



((710, 0, CT30, CT40, 050) € I? 



6 , 4 (7ioCr| 

— (730040 ~ 01O0'5O + - 



5 5 1-crf 



+ AcrioY^l -0? o = °J 
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These regular equilibria, in the original system, correspond to periodic motions of the disk in 
which the point of contact describes a circle and the center of mass stands at constant height. 
These motions are contained in the set of steady motions of the rolling disk, according to 
Routh's terminology [33, 35]. They have received an extensive treatment in [18], although by 
using the system (26). 

Now, both of the systems (23) and (26) admit two first integrals related to the solutions 
(in the sense explained in Section 3.1) of the non-autonomous linear system 

da 3 „ da 4 2 cr 3 

- — = -2er 4 , - — =--■ ( 2S ) 

do\ do\ 6 1 — a\ 

which can be written in matrix form as 

d f "\-( a -:-)(? i. (29) 



d<7l V CT 4 / V ~3W? / V ^ 

This equation is the same as Eq. (69) of [18], where its solutions have been studied in great 
detail, including their asymptotic behaviour. 

However, the important point for us is that the systems (23) and (26) are good candidates 
to be formulated as Hamiltonian systems with respect to Poisson structures of the type 
described in Theorems 2 and 1, respectively. Let 9 — d(f> and 6\, #2 be the 1 -forms, defined 
in M (resp. V) by 

2 cr 3 

Applying the results of Sections 3.1 and 3.2 to these 1 -forms, we have 



h = 2(T4 do\ + das , 82 = — j do~i + da^ . 



Proposition 2 The bivectors of the form A = — A i2 U A V, defined in M, where 

d rr d n d 2 cr 3 d 

U =7^> V =7^- 2(T ^--iT^7^> (30) 
0<T2 oo\ a<73 3 1 — af 004 

and A12 € C°°(M) is a non-vanishing function, are Poisson tensors of rank two in M. 

The vector field X in M, whose integral curves are the solutions of (26), is a Hamiltonian 
vector field with respect to the Poisson bivector A with the specific function A12 = 1 — a\ and 

Hamiltonian function E given by (27), i.e., X = A (dE) in M. 

Proposition 3 The bivectors A = f[(Z(j>)U A V + Y A Z], defined in V C K 5 , where U and 
V are given by (30), Z = d/da^, Y = (U cf))V — (V <j))U , and f € C°°(2?) is a non-vanishing 
function, are Poisson tensors of rank two in V, except in the set of singular equilibria, where 
they vanish. 

The vector field X in T>, whose integral curves are the solutions of (23), is a Hamiltonian 
vector field with respect to the Poisson bivector A with the specific function f = 1 and 
Hamiltonian function E given by (24), i.e., X = A$(dE) in V. 

Both Propositions can be proved by direct computations. 

The Poisson Hamiltonian structure of the systems (23) and (26) could be used to have 
an interpretation of their geometry. For example, the invariant submanifolds mentioned in the 
analysis of the reduced vector field (26) in [18], could be understood as the symplectic leaves 
of the rank-two Poisson structure(s) A of Proposition 2. 
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4.2. Routh 's sphere 

For this example we will follow the treatment and use some of the results of [17], see details 
therein. This problem has been treated as well, e.g., in [4, 7, 20, 33, 35]. Consider a sphere of 
mass m and radius r with its center of mass at a distance a (0 < a < r) from its geometric 
center. The line joining both centers is a principal axis of inertia, with associated moment of 
inertia I 3 . Any axis orthogonal to the previous, passing though the geometric center, has an 
associated moment of inertia I± . This sphere is supposed to roll on a horizontal plane under the 
influence of a vertical gravitational field of strenght g. The resulting non-holonomic system 
has as well two symmetry groups. One is again the group E(2) consisting of translations 
in the horizontal plane and rotations about the vertical axis. The other is the S 1 symmetry 
consisting of rotations about the principal axis of inertia which joins the center of mass and 
the geometric center of the ball. 

Again, after these symmetries have been reduced out by a similar procedure to that of 
the rolling disk, it is obtained a system giving the evolution on the reduced orbit space, which 
coincides with that of the rolling disk. 

Therefore, let us consider a reference ball as the one described, with the geometric 
center at the origin, and the center of mass at the point — ae 3 , where e 3 denotes the vertical 
unitary vector in this fixed frame. The position and attitude of the moving ball is given by 
transforming the position of the reference ball by means of a translation a (e.g., of the center 
of mass) and a rotation A. We denote by s the vector in the fixed sphere such that rotated 
by A gives the vector in the moving sphere pointing from the center of mass to the point of 
contact. The unitary vector u in the fixed frame is the pre-image of — e-$ under the rotation A, 
u = —A~ 1 es. The relation between u and s is 03 = (s, u). The components of the angular 
velocity uj of the ball will be denoted by (io\ , u>2 , W3). 

Following [17], after the reduction of the mentioned E{2) symmetry, the equations of 
motion read 



The second of Eqs. (31) expresses again the non-holonomic constraint of rolling without 
slipping. 

Now, the reduction of the S 1 symmetry is performed in an analogous way as in the case of 
the rolling disk, see Section 4.1, where (ui, u 2 , u 3 ) denote as well the components of u. The 
S 1 action consists of rotating both vectors u, lo simultaneously, with respect to the principal 
axis joining the geometric and mass centers. This action is not free, since S 1 leaves invariant 
pairs of points of the form ((0, 0, ±1), (0, 0, LO3)). The corresponding set of invariants is again 
(19) with the relations (20). Thus, the reduced orbit space M is the semialgebraic variety of 
IR 5 described in the previous example of the rolling disk, with the same notations. 

However, the reduced system reads now, using (31), 



— {Ilo + ms x (oj x s)) =/wxw + m— x (lo x s) 




(31) 



(32) 



<7l = (72 

T(<7i)<72 = (/3 + mr 2 + mraai) 0-3(74 — mga(l — a\) 

— a$(mra + (h + ma 2 + mr 2 )a\ + mraa 2 ) 
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o 3 = - h ^ g4 , (7 3 + mr 2 + mraoi) (33) 

O2O4 , , . . 

o- 4 = - mr- — -(7 3 a + r(7 3 - hpi) 

T(ai)&5 = — 2mra<72<J5 — 2mgaa 2 

2 7 3 + mr 2 + mraoi 
- 2mr (7 3 - 7i) — — r cr 2 (T 3 cr4 , 

where P(oi) = 7i7 3 + mr 2 7i(l - of) + m7 3 (a + roi) 2 and T(oi) = 7i + rw 2 + ma 2 + 
2mraoi . The reduced energy is 

75 = i(T(oi)o 5 + (7 3 + mr 2 )o 2 - mr^s + cn^) 2 ) + ma(gu x - ra 3 a 4 ) ■ (34) 

These are eqs. (23), (24) and (25) in [17]. In this case the expressions (33) and (34) make 
sense for all T> = R 5 . We will consider this extended domain for the vector field X whose 
integral curves are given by (33) and for the reduced energy function E. 

Restricting ourselves to points in M with |oi | < 1, we find again the smooth submanifold 
M C M given by (25). The integral curves of the projected vector field X are the solutions 
of the system ( [17], eq. (38)) 

Oi = 02 

T(ai)&2 = (7 3 + mr 2 + mraa\)a 3 a4, — mga(l — of) 

'mra + (h + ma + mr )o\ + mraaf) (35) 



T u 2 u 4 2 

& 4 = - mr J 2 ,' 74 , (ha + r(7 3 - 7i)oi) 

and the restricted reduced energy E is 

\ ( T ( gl ) yj~3r + ( J 3 + ™r 2 )o| - mr 2 (o 3 + 0!0 4 ) 2 ^ + ma(<?Oi - ro 3 o 4 ) . (36) 

The reduced vector field X satisfies X(E) = and X(<£) = in D, and = in 

M. Moreover, X has a family of singular equilibrium points, belonging to the singular set II± , 
given by {(±1, 0, 0, o 4 , 0) | o 4 e M}, which physically correspond to the spinning of the ball 
about its symmetry axis when it is vertical (then the reduced energy becomes ^7 3 o| ± mga). 
It has as well a family of regular equilibria given by the set of constants 

{(oio, 0, o 30 , o 4 o,o 5 o) G M 5 I &(oio, o 30 , o 40 , o 50 ) = 0} , 

where 6(01, o 3 , 04, 05) = (7 3 + mr 2 + mra.o\)o 3 a\ — mga(l — a\) — <r 5 (mra + (7i + 
ma 2 + mr 2 )ai + mraa 2 ). These regular equilibria, in the original system, correspond to 
periodic motions of the ball in which the point of contact describes a circle and the center of 
mass stands at constant height. 

In this case, both of the systems (33) and (35) admit two first integrals related to the 
solutions (in the sense of Section 3.1) of the non-autonomous linear system 

c?o 3 7 3 (7 3 + mr(r + aoi))o 4 do 4 mr(l\ro\ — 7 3 (a + roi))o 4 



rfoi P(oi) ' da x P(oi) 



(37) 
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Thus, the mentioned systems are other good candidates on which to apply the Poisson 
approach of Section 3. Let 9 = d(f> and 9\, 9 2 be the 1-forms, defined in M (resp. V) 
by 



n = 



->2 = — 



mr(I\ra\ — I 3 (a + ro\))o4 



dai + da 3 , 



da\ + do 4 . 



We have the following results: 

Proposition 4 The bivectors of the form A = — A i2 U A V, defined in M, where U = 8/802, 

v _ _d I 3 (I 3 + mr(r + aox^Oj d | mrjhroi - I 3 (a + ro 1 ))<74 d ^ 

doi P{ a i) 9o 3 P{°~i) do4 ' 

and A12 € C°°(M) is a non-vanishing function, are Poisson tensors of rank two in M. 

The vector field X in M, whose integral curves are the solutions of( 35), is a Hamiltonian 
vector field with respect to the Poisson bivector A with the specific function A12 = (1 — 

ol)/T(oi) and Hamiltonian function E given by (36), i.e., X = A\dE) in M. 



Proposition 5 The bivectors A = f[(Z<f>)U A V + Y A Z], defined in V = R 5 , where U and 
V are given as in Proposition 4, Z = d/do 5 , Y = {UcfyV - (V<j))U, and f e C°°(V) is 
a non-vanishing function, are Poisson tensors of rank two in V, except in the set of singular 
equilibria, where they vanish. 

The vector field X in T>, whose integral curves are the solutions of( 33 ), is a Hamiltonian 
vector field with respect to the Poisson bivector A with the specific function f = 1 /T(o~i) and 
Hamiltonian function E given by (34), i.e., X = A$(dE) in V. 

Both Propositions can be proved as well by direct computations. 

In this case, the equations (37) can be explicitly integrated in an easy way. From the 
second of these equations we have the relation er 4 y/P(a\) = k. Substituting into the first, 
we can also integrate to obtain the relation I\ro 3 + I 3 (a + ro\)o4 = j. The constants 
k,j are integration constants (essentially the initial conditions of the system (37)). These 
two expressions are the desired first integrals (Casimir functions of the preceding Poisson 
structures). The second of them is known as Jellet's integral, see [17,20] and references 
therein, see also p. 184 of [7]. 

The invariant submanifolds thoroughly studied in [17], could be interpreted in this 
framework as the symplectic leaves of the rank-two Poisson structure(s) A of Proposition 4, 
determined by the level sets of the first integrals j and k. 



4.3. Ball rolling on a surface of revolution 

For this example we will follow the treatment and use some results of [25], see therein for 
more details. This problem has been treated as well, e.g., in [8, 33, 35]. In particular Routh, 
in the last of these references, noticed the existence of two integrals of motion given by a 
system of two linear differential equations, solved them in special cases, and described a 
family of stationary periodic motions together with a necessary condition for their stability. 
Later, in [42], it has been shown that the condition is also sufficient. Both of [25] and [42] 
prove that the corresponding reduced system has integral curves consisting of either periodic 
orbits or equilibrium points. 
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Consider a homogeneous ball of mass m, radius r and moment of inertia M with respect 
to any principal axis. The ball rolls without slipping on a surface of revolution, under the 
influence of a vertical gravitational field of strenght g. We take the origin of coordinates at 
a point of the axis of symmetry of the surface (the intersection of this axis with the surface 
at its vertex), and we consider a horizontal plane passing through it. We parametrize the 
position of the center of mass of the ball by its coordinates (01, a?) on this horizontal plane, 
and its height will be parametrized via the smooth profile function (p : M M of the surface, 

a 3 = (p (y\fa{ + a£j. Note that not all surfaces of revolution can be parametrized well 
in this way, e.g., the cylinder, which requires a separate treatment, see Section 4.4 below. 
We will assume that tp is a smooth even function, thus we will have that ^ 2fc+1 )(0) = 0, 
fc = 0, 1, 2 . . . We denote by (ui,u>2, 0J3) the components of the angular velocity vector lo of 
the ball, and (71, 72, 73) the components of a unit vector 7 normal to the surface at the point 
of contact (directed towards the center of the ball). The unit vector in the vertical coordinate 
axis is es. 

The equations of motion can be easily computed by the classical equations of the 
variation of the angular momentum, and implementing the non-holonomic constraint of non- 
slipping of the point of contact, i.e., that its instantaneous velocity vanishes. They read (with 
respect to the center of mass of the ball, compare with eqs. (5), (7) of [25] and Section 2 
of [8]) 

M^- — mr 2 ( x 7) J x 7 — mgre 3 x 7 = , 

CLL \ LLL J 

a - r(uj x 7) = . (39) 
The total mechanical energy of this system is then 

H = i((M + mr 2 )(uj ■ lo) - mr 2 (j ■ uj) 2 ) + mga 3 , (40) 

and is a first integral for the system (39). 

The system (39) and the energy (40) admit a further reduction of the S 1 symmetry 
consisting of rotations of the system about the vertical axis, and thus rotating both of lo 
and 7 simultaneously. This action, as in the previous cases, is not free, since the isotropy 
subgroup of pairs ((0, 0, 1), (0, 0, 0J3)) is S* 1 (these pairs correspond to motions of the ball 
spinning around the vertical axis when being at the vertex of the surface), and we will use 
again invariant theory in order to perform the reduction, but now as it has been done in [25]. 
First of all, we define the vector v and the scalar w as follows: v = r(uj x 7), w = —r(uj ■ 7). 
Then, a full set of invariant polynomials, which parametrize the orbit space of the S 1 action, 
is 

Pi = \( a l + a i) 1 P2 = CLlV! + a 2 v 2 , p 3 = aiV2 ~ «2«1 , 

p 4 = w, Pn = \{vl + vl), (41) 



2 



with the relations 



pl+pl~4piP5 = 0, pi>0, p 5 >0. (42) 
Therefore, the completely reduced orbit space P is now the semialgebraic variety of M 5 

P = {( Pl , ...,ft)6l 5 I 0(p) = 0, Pl > 0, P5 > 0} , (43) 

where now0 e C°°(M 5 ) is the polynomial function 4>(p) = p\+p\— ^PiPb- P is not a smooth 
submanifold of M 5 , because the previous S 1 action is not free. Instead, P is homeomorphic 
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to a cone in R 4 times R [25], which can be easily seen from the relation (f>(p) = when it 
is written as p\ + p\ + (p 1 — p 5 ) 2 = (pi + p$) 2 . The vertex of the cone is determined by 
P2 = P3 = Pi = P5 = 0, therefore the singular points of P are 



n = {(o,o,o, P4 ,o) e p | p 4 e R}. 

Calculating the time derivatives of the invariants, using (39), and the relations 



7i 



ai 



/ Wl^/l + if' 2 



72 



02 



73 



1 



(we will use the notation ip = <p(y/2pi), <p' = <p'(y/2pi) and <p" 
follows) we arrive to the system in the reduced orbit space P 



ip"(y/2pi) in what 



Pi = P2 



1 



P2 



P3 



Pa 



P5 



1 + ip' 
M 

ar z 1 



M 

OLT 2 



TP3P4- 



/ 2pi(p' + 2p 5 



Wi 



P2P3 
' 2pi 

P2 



(44) 



1 



/2fi 



and the reduced energy 

E 



f 1 (M 

' (_ 2pi \ otr z 



Pi 



iWi 



iWi 



tng v 

/2p~! 



a 



2P5^- 

2pi 



M 2 atp'- 
_p 4 + ap5 + _ 



■p 2 + mg(p , 



(45) 



where a = M+ f r . These are the equations found in Lemmas 2.2 and 2.3 (i) of [25]. 

We observe that the right hand sides of (44) and (45) make sense in an open set V of R 5 
larger than P, namely T> = R 5 \{(pi,p2,P3,P4,P5) | Pi < 0}. This is due to the fact that they 
are defined in the limit p\ — ► + , because of the assumption that the odd-order derivatives at 
of ip vanish. (For points strictly in P with p\ = this assumption would not be necessary, 
since these points also have p 2 = 0, p^ = 0). We will consider the enlarged domain T> for the 
vector field X whose integral curves are the solutions of (44), and also for the reduced energy 
(45), compare with p. 500 of [25]. 

The regular stratum of P, i.e., P\H, can be covered by two charts [23], whose 
corresponding neighbourhoods can be chosen to be the smooth open dense subsets Pi , P 2 C 
P given by 



Pi = \ (P!,P2,P3,Pi,P5) e 



{pi,P2,P3,Pi,P5) e 



P5 



Pi 



Pi 



■Pi 



4pi 



Pi > o 



(46) 



P\ 



P% 



4p 5 



Pi 



> 



However, for our purposes here, it will be enough to consider just Pi, in order to 
endow it with Poisson structures of the type described in Theorem 1, which afterwards could 
be compared with the Poisson structure given originally in [8]. The procedure for P2 is 
analogous. Thus, we will consider the induced vector field X and energy E on Pi, which can 
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be found from (44) and (45) by solving for p 5 . The integral curves of X are the solutions of 
the system 

Pi = P2 

1 f M Lp' mg r— , pl+pl 2 f' ( „ <p' 



1 + Lp' [ ar 2 \/2~Pi a 2pi \/2pi V V2p 



M 

= — 2P2P4— (47) 



P4 = 



ar" 1 + Lp'- 

P2PZ ( f" <fi' 



2pi VI + ^' 2 V2w 
meanwhile 



E =2^ pl+ at %r + + mg * • (48) 

Now, the reduced vector field X satisfies X(E) = and X(^) = in V, and = 

in Pi. The vector field X has a family of singular equilibrium points consisting of the 
singular set II, that is, {(0, 0, 0,p 4 , 0) | p 4 e R}, which as already mentioned, correspond 
to the spinning of the ball about the vertical when being at the vertex of the surface (then the 
reduced energy becomes j^p\ + mgtp(0)). X has as well a family of regular equilibria given 
by the set of constants 

|(Pio,0,p 3 o,P40,P5o) € V 2p 50 - ^-^2p w ip'(^2pi ) - ~~2 ^ [ffi^ P3oP40 = °} 

These regular equilibria correspond in the original system to rotations of the ball along a 
parallel of the surface of revolution at constant height. 

In addition, both of the systems (44) and (47) admit two first integrals of motion related 
to the solutions (in the sense explained in Section 3.1) of the non-autonomous linear system 
(see also [25], Lemma 2.3. (ii)) 

dps M Lp" d P4 P3 ( <p" V 1 \ (aq\ 



dpi ar 2 1 + Lp' 2 dpi 2pi \ 1 + ip' 2 \/2pi j 
Let #o = d(f> and 6i, 82, be the 1 -forms, defined in Pi (resp. 2?) by 

M Lp" p 3 ( Lp" Lp' \ 

Oi = — tP4— — j^dpi - dp 3 , 6 2 = — — — H -== dpi + dp 4 ■ 

ar 2 1 + ip' 2 2pi \1 + tp' 2 \/2pi J 

We have the following results, applying the Theorems of Section 3, which can be proved by 
direct computations: 

Proposition 6 The bivectors of the form A = — A12 U A V, defined in Pi, where 

U=— V= 9 I M p, ^ 9 P3 ( V " fL)l (50) 

dp 2 ' dpi ar 2 1 + tp' 2 dp 3 2pi \ 1 + tp' 2 \/2p x ) dp^ ' 

and A12 € C°°(Pi) is a non-vanishing function, are Poisson tensors of rank two in Pi. 

The vector field X in Pi, whose integral curves are the solutions of (47), is a 
Hamiltonian vector field with respect to the Poisson bivector A with the specific function 
A12 = 2pi/a(l + Lp' 2 ) and Hamiltonian function E given by (48), i.e., X = k\dE) in Pi. 
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Proposition 7 The bivectors A = f[(Z<j>)U A V + Y A Z], defined in V C R 5 , where U and 
V are given by (50), Z — d/dp^, Y = (U (f)V — (V 4>)U , and f G C°°(T>) is a non-vanishing 
function, are Poisson tensors of rank two in T>, except in the set of singular equilibria, where 
they vanish. 

The vector field X in T>, whose integral curves are the solutions of (44), is a 
Hamiltonian vector field with respect to the Poisson bivector A with the specific function 
f = l/2a(l + Lp' 2 ) and Hamiltonian function E given by (45), i.e., X = A*(dE) in V. 

Remarks For the present case, a Poisson structure analogous to one of the structures A of 
Proposition 6 has been found, to the best of our knowledge, by the first time in [8], see 
their equation (3.11) for A = 0. In fact, up to a rescaling, they are the same, by using the 
identifications 



X 2 = —p i y / 2p 1 , x 3 = — = , Xi = ar == , 

r if' r ^/2^V 1 + y/2 yf2p{ 

Xl = 7^' f{xi) = ✓ ' 

The (local) Poisson bivector found in [8] for this case reads in their coordinates 

(xi,X2,x 3 ,Xi) as 

2 / 8 f'( Xl ) d \ 2 X 2 d \ d 



dx\ x\ dx2 ) fi x i) dx 3 J dx 



''4 



which, in particular, is also of the type described in Example 1 . Therefore, multiples of this 
bivector are again Poisson bivectors and hence, the rescaling introduced in [8], by means of 
an invariant measure, in order to render the reduced system Hamiltonian, is unnecessary. 

On the other hand, Hermans in [25] has not noticed the existence of any of these Poisson 
structures of rank two but he constructed a closed 2-form, with domain contained in Pi, which 
vanish in a set containing the set of regular equilibria, but has rank four otherwise. For this 
construction, which uses non-holonomic reduction [5], it is indeed necessary to rescale the 
original reduced vector field, see Section 4. 1 of [25]. 



4.4. Ball rolling on the interior of a cylinder 

In this Section we will treat the special case of a ball rolling inside of a cylinder, which cannot 
be parametrized as in Section 4.3. In contrast with the general case, this case is completely 
and explicitly solvable, as it is well known, see, e.g., [4, 8, 32, 33]. However, we will give an 
independent treatment. 

For this specific system, we will easily find a family of Poisson structures of rank two, 
generated by two of them, with respect to which the reduced system is Hamiltonian with the 
reduced energy as Hamiltonian function. 

Consider therefore the case of the ball rolling inside a surface of revolution, with the 
following variations: the center of mass of the ball will be parametrized by the vector a, with 
cylindrical coordinates (p cos 9, p sin 6, z), where p is the radius of the cylinder on which the 
center of mass of the ball moves, and z is the height with respect to the gravitational energy 
reference point. The normal vector 7 reads then 7 = — (cos#, sin#, 0). The system (39) 
becomes in the coordinates (w\, ui2, W3) and (0, z) 

m ( g r \ 
Co\ = — — I — oj 3 (uji cos 8 + 0J2 sin 9) sin 9 , 
a \r p J 
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d>2 = ( — I — ^3(^1 cos 9 + 0J2 sin#) ) cos 8 , (51) 

a \r p J 

r 

w 3 = , 9 = — 0J3 , z — r(u2 cos 9 — uji sin 9) . 
P 

Likewise, the energy (40) reads now 

H = i{(M + mr 2 )(uj ■ to) - mr 2 (Lj 1 cos 9 + uj 2 sin6>) 2 } + ragz , (52) 

which is conserved by the system (51). Obviously, uj 3 is a first integral of the system as well. 

Let us consider now the system obtained after the reduction of the S 1 symmetry of 
rotations of the whole system about the vertical axis, as in the general case. Although now 
the S 1 action is free, we will use again invariant theory in order to perform the reduction. 
Consider the invariants similar (but not equal) to (19): 

<7i = z , cr 2 = 71^2 — 72^1 = —^2 cos 9 + uj\ sin 9 , 

C73 = 71 wi + 72^2 = — wi cos 9 — ul>2 sin 9 , 04 — uj 3 , (53) 

which in this case can be regarded as coordinates on R 4 . Then, the reduced system for 
(<7i, cr 2 , 03, 04) reads 

Ma '4 rag . r 

ai = -ra 2 , a 2 = a 3 H , a 3 = — a±a 2 , 04 = , (54) 

arp ar p 

which preserves the reduced energy 

E = ^{mr 2 {al + a\) + M (a\ + a\ + a\)} + mga 1 . (55) 

The reduced vector field X in the reduced space reads then 

d ( Ma t . mg\ d r d 
X = ~ra 2 -^ h a 3 H 04,0?.— — , (56) 



' dai \ arp ar J 802 p da 3 

and we have X(E) — in all points of M 4 . The general solution of (54) can be given 
explicitly. It reads 



a 1 (t) 


= *i(0) 


*2(t) 


= <7 2 (0) 


<*»(*) 


= _f£ 


V2 




= 04 , 



r 



-{tr 2 (0)(l - cosy/vi^t) + \/v\Vi 0-2(0) sva^/v^t} 

Vll>2 

„ m , *2(0) 



sin yjviv 2 t (57) 



~ ^ cos yJv\W2.t — \ h^- a 2 (0) sin yJv\Vit 



V2 \ v 2 



where we have defined the constants v\ — ra 4 /p, v 2 = Mv\jar 2 and a g — mg/ar. It is 
clear that the reduced system, if 04 7^ 0, has integral curves consisting of either periodic orbits 
or equilibrium points, belonging to the set {(010, 0, — mgp/Ma4o, 040) G K 4 | 010 ; 040 € 
K, 040 7^ 0}. These equilibrium points correspond to rotations of the ball inside the cylinder 
at constant height, as in the general case. On this occasion, the reduced system can be 
reconstructed easily to the complete system, thus the general solution of (51) is 

= 0- 2 (t)sin0(i) - a 3 (t) cos 9(t) , 
lo 2 (t) = -a 2 (t) cos 9(t) - 03 (t) sin 6(t) , (58) 
w 3 (i) = ct 4 = w 3 (0) , z(t) = a 1 (t), 
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where 0(t) = 9 — v\t and a^t), i = 1, 2, 3 are given by (57). If we denote Wi(0) = u> w , 
w 2 (0) — uj 2 o, we have the relations for the initial conditions 

cr 2 (0) = ujio sin 6»o - lu 2 q cos0 o , ^(O) = <*g — ^2(^10 cos# + ^20 sin# ) • 

The complete system, when ^3 7^ 0, is then isochronus with two frequencies, the motions 
being periodic (relative equilibria, projecting to equilibrium points in the reduced space) or 
quasi-periodic, otherwise. The solutions with ui 3 = correspond to falling motions of the ball, 
rolling along a vertical generatrix of the cylinder. The explicit expression of these solutions is 

wi(t) = u>io + ta g sin a , lu 2 (t) = w 20 - t a g cos O , u 3 (t) = , 
6{t) = O , z(t) = z - 7} r<J g t2 + rt ( w 20 cos O - oj 10 sin O ) . 

We will treat now the question of writting the vector field X, given by (56), as a 
Hamiltonian vector field with respect to a Poisson structure of rank two, with Hamiltonian 
function E. 

We first observe that the reduced vector field X is annihilated by the 1 -forms 

9i = da 4 , 2 = -da\ + da 3 , 3 = ar 2 a2a 4 da2 + (Ma 3 a 4 + mgp)da 3 , 

and then, it is easy to apply Theorem 1, to obtain the following results: 

Propositions The bivector Ai = ^ A ±X = a (r-^ + ^af;) « a Poisson 

bivector on K 4 of rank two such that A*(0i) = A* (# 2 ) = 0. Likewise, The bivector A 2 = 



JL a AX = — —fr-A \( ^a 3 + ma)d_ £(74(72 

mg ao\ mg oai \ arp ° ar j oa? p * * 



da 3 



is a Poisson bivector on } 



of rank two such that A 2 (#i) = A 2 (0 3 ) = 0. In addition, we have X = A\(dE) = A\(dE), 
where X is given by (56) and E by (55). 

Now, the Pfaff systems "6>i = 0, 2 = 0" and "#1 = 0, 3 = 0" can be easily integrated, 
giving non-trivial Casimir functions of Ai, A 2 , and first integrals of X: 

Proposition 9 We have kerA* = spanjrfci, c£c 2 }, and kerA 2 = spanjrfci, dc 3 }, where 
ci = a 4 , c 2 =a 3 - Zfay and c 3 = r -fa 2 + (^a 3 + ^) a 3 . 

As a consequence, we have that the reduced vector field X has in principle four first 
integrals, namely, E, a, c 2 and c 3 , but clearly, they form a functionally dependent set. 
However, for example, we have that {E,C2,c 3 } is generically an independent set of first 
integrals, although in the equilibrium points one becomes dependent of the other two. In the 
falling motions, 04 = 0, therefore Ai and A 2 become proportional. 

Incidentally, we also observe that A{(dc 3 ) = -^X and A 2 (rfc 2 ) = --^X. In 
addition, the Poisson bivectors Ai, A 2 are compatible in the sense that their Schouten- 
Nijenhuis bracket vanishes, [Ai, A 2 ] = 0, which can be checked, e.g., using the properties 
(1). Thus we have the following result: 

Proposition 10 The pencil of bivectors A\ = (1 — A)Ai + AA 2 consists of Poisson bivectors 
of rank two such that X = A^ x (dE) for all A G M. Moreover, the functions c\, c 2 a = 
(1 — X)E — ar 2 c 3 and c 3 \ — Xra^E/p + mgr c 2 are (functionally dependent) Casimir 
functions of A\ (and therefore, first integrals of X) for all A G K. 

Proof 

That the rank of A\, for all A € 1, is two, is obvious when one realizes that it does not 
contain terms on an d therefore the rank must be an even number between zero and four. 
The other statements are a matter of computation using the above observations. 
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5. Conclusions and outlook 

We have shown the form of certain Poisson structures of rank two with respect to which 
certain reduced problems of non-holonomic mechanics become Hamiltonian. We have shown 
that in K 4 and R 5 , from an algebraic point of view, these Poisson structures are defined, up 
to a factor function, by the choice of the kernel of bivectors on these spaces to be generated 
by 1 -forms of a specific type. Such 1 -forms define integrable codistributions in the sense of 
Frobenius, and are chosen in order to accommodate and generalize the systems of first order 
non-autonomous differential equations which appear after reduction in certain non-holonomic 
mechanical systems, whose solutions are related to first integrals of such reduced systems. 

We have applied the theory to the cases of the rolling disk, the Routh's sphere, and 
the ball rolling on a surface of revolution, explicitly recovering as a particular case some 
results of [8]. Thus, we have shown that the framework suggested by Borisov, Mamaev and 
Kilin [7, 8] can be improved along the lines discussed, namely, that those reduced systems 
need no rescaling to become Hamiltonian with respect to a Poisson structure of rank two, and 
that the domain of definition of the Poisson structures introduced therein can be extended, 
including even the set of singular equilibria of the reduced systems. A natural question is 
whether a similar approach could be used in other non-holonomic systems, maybe of higher 
dimension. 

However, there are more fundamental points still to be better understood. For example, to 
what extent the mentioned Poisson structures can be useful to investigate the intimate nature 
of these and maybe other non-holonomic systems, for example in order to characterize their 
integrability properties [3,21,22]; see also the recent work [23]. Another question could be to 
clarify the origin of the system of differential equations giving the conservation laws for the 
mentioned reduced non-holonomic systems, see also [4, 10, 19, 32,40] and references therein. 
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